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Abstract 



We compute the quadrupole mode of the gravitational waveform of inspiralling compact binaries 
at the third and half post-Newtonian (3.5PN) approximation of general relativity. The computation 
is performed using the multipolar post-Newtonian formalism, and restricted to binaries without 
spins moving on quasi-circular orbits. The new inputs mainly include the 3.5PN terms in the 
mass quadrupole moment of the source, and the control of required subdominant corrections to 
the contributions of hereditary integrals (tails and non-linear memory effect). The result is given in 
the form of the quadrupolar mode (2, 2) in a spin- weighted spherical harmonic decomposition of the 
waveform, and may be used for comparison with the counterpart quantity computed in numerical 
relativity. It is a step towards the computation of the full 3.5PN waveform, whose knowledge is 
expected to reduce the errors on the location parameters of the source. 
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I. INTRODUCTION 



The first ever direct detection of gravitational waves might occur around 2016, when the 
current network of ground-based laser-interferometric detectors (such as VIRGO and LIGO) 
will have been upgraded to a higher sensitivity. Inspiralling and merging binary systems 
composed of black holes and/or neutron stars are among the most promising sources for 
those detectors. Though the gravitational waves are extremely weak, the large number of 
predictable cycles in the detector's bandwidth will enable the on-line detection and later 
the parameter estimation of the signal. The detector output will be cross-correlated with a 
number of copies of the theoretically predicted signal or template, corresponding to different 
signal parameters. The need of a faithful template bank has driven the development over the 
last twenty years of both accurate approximation methods and powerful numerical schemes 
in general relativity. For a review of gravitational wave detectors and future gravitational 
wave astronomy, see pQ. 

The templates for the inspiral and merger of two compact objects (say black holes) are 
computed by combining a high order post-Newtonian approximation for the early inspiral 
phase [2], with a full-fledged numerical integration of the field equations for the late inspiral 
and ringdown phases [3H7j . The post-Newtonian and numerically-generated results are then 
matched together with high precision, yielding the full gravitational waveform, including all 
amplitude and phase modulations, described either analytically and/or numerically [8T4T5]. 

Previous post-Newtonian works [ToT - [2T)] have provided the waveform including all its har- 
monics (i.e. beyond the dominant harmonic at twice the orbital frequency) up to the 3PN 
orderQ In applications to data analysis the full waveform should be used for detection and 
parameter estimation up to the maximum available post-Newtonian order. In particular, it 
should further improve the angular resolution and the distance measurement of the system 
for massive enough binaries [211 E2]. Now, the quadrupole mode at twice the orbital fre- 
quency, having (£, m) = (2, 2) in a spin-weighted spherical harmonic decomposition, is the 
dominant one in the sense that it gives the only time- varying contribution to the waveform 
at the dominant orderj^] It is also the one which is computed with the best precision in 
most numerical simulations. In the present paper we extend the previous works [T6TI2T)] by 
computing the dominant (2, 2) mode at the next 3.5PN order. Our result agrees with the 
one derived within the black-hole perturbation theory in the limit where the binary mass 
ratio goes to zero [231 El]- The completion of the other modes at the 3.5PN order does not 
appear to be straightforward, notably regarding the (2, 1) mode, and will be left for future 
investigation. 

The computational basis is the multipolar post-Newtonian wave-generation formalism, 
which has two different aspects. First, it constitutes a general method applicable to extended 
sources with compact support, which combines a mixed post-Minkowskian and multipolar 
expansion for the field outside the source [25-28J, with a post-Newtonian expansion for the 
field inside the source [29|. Second, it addresses the problem of applying this method to 
point-particles modelling compact objects [301 IS], in which case it crucially requires a self- 
field regularization |32j . In the present article we shall mainly focus on the results and refer 
to previous papers for full details on this formalism (see also [20] for a summary)]^] 

1 As usual we refer to nPN as the post-Newtonian terms with formal order 0(c~ 2n ) relative to the Newtonian 
acceleration in the equations of motion, or to the quadrupole-moment formalism for the radiation field. 

2 The (£, m) — (2, 0) mode contributes at the dominant Newtonian order, but only in the form of a non- 
oscillating (DC) term. 

3 In the following we shall refer to Ref. [20] as Paper I. 
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Our plan will be as follows. In Sect. [IT] we recall the equations of motion of non-spinning 
compact binaries on quasi-circular orbits up to 3.5PN order. In Sect. |III| we remind some 
basic definitions for computing the gravitational waveform and associated polarization modes 
for planar compact binaries. In Sect. |IV| we express the dominant radiative quadrupole 
moment in terms of the source multipole moments of a general isolated source up to 3.5PN 
order. In Sect. |V]we give the expressions of those source moments fully reduced in the case 
of circular compact binaries at the same accuracy level. Finally the result for the (2, 2) mode 



at 3.5PN order is presented in Sect. VI The source quadrupole moment at 3.5PN order for 



non-circular binaries in a center-of-mass frame is relegated to Appendix |Xj 



II. QUASI-CIRCULAR BINARY AT 3.5PN ORDER 

An inspiralling compact binary of non-spinning compact bodies is modelled as a system 
of two particles solely described by their masses mi and m 2 . The orbital plane is spanned 
by the relative position of the particles x = yi — y 2 and the relative ordinary velocity 
v = dx/dt. The unit vector normal to this plane is given by £ = x x v/\x x v\ (we 
assume a non-radial orbit) and is constant in the absence of spin effects. Introducing the 
unit separation direction n = x/r, where r = \x\ is the separation distance, and posing 
A = t x n we have 

x = rn , (2.1a) 
v = rn + rco\, (2-lb) 
dv 

— — — (r — r u 2 ) n + (r uj + 2r uS) A , (2.1c) 

which defines the orbital frequency u related in the usual way to the orbital phase <p by 
u — (p. The motion of the binary follows a quasi-circular orbit decaying by the effect of 
radiation reaction starting at 2.5PN order. Using the facts that r = 0(c~ 5 ) and u = 0(c~ 5 ), 
and noticing that f = 0(c~ 10 ) is of the order of the square of radiation reaction effects, we 
obtain the quasi-circular acceleration as 

dv o (Cj r\ _ / 1 \ , 

-co 2 x+ (- + 2-)v +0 — . (2.2) 



dt \lu r/ \c 10 1 

The conservative part of the dynamics is given for circular orbits by the expression of 
the orbital frequency u in terms of the binary's separation r up to 3PN order. This result 
has been obtained in harmonic coordinates [33H39] and in Arnowitt-Deser-Misner (ADM) 
coordinates [I01 - H2] . For the present work r is the orbital separation in harmonic coordinates, 
and from it we define the post-Newtonian parameter 

Gm 

7 = — • 2.3 

rcr 

Our mass parameters will be the total mass m = nix + m 2 and the symmetric mass ratio 
v = m 1 m 2 /m 2 . The orbital frequency is then given by the 3PN "Kepler's law" 

c 2 = ^ji + 7 (-3 + z/)+ 7 2 (e + ^ + ^ 2 



3 



+ 7^ 



-10 



75707 41 



840 



64 



-7T 



22 In 



o/ J 



19 5 

z/ H v 1 - 

2 



(2.4) 



We neglect the 4PN and higher terms (9(7 4 ). The logarithm at 3PN order comes from a 
Hadamard self-field regularization scheme [331 E3] and involves a regularization constant r' 
specific to harmonic coordinates. This constant r' will disappear from our physical results 
in the end. By inverting (2.4) we obtain 7 in terms of the alternative parameter 



x 
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(2.5) 



which is an invariant in a large class of coordinate systems^] including for instance the 
harmonic and ADM coordinates. At 3PN order we have 
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(2.6) 



The dissipative radiation reaction part of the equations of motion can be computed by 
balancing the change in the orbital energy E with the total energy flux J 7 radiated by the 
gravitational waves, dE/dt = — J 7 . Up to 3.5PN order, which means using the orbital energy 
and energy flux at 1PN relative order, this gives 



64 Gm 



r — — - 



v 7 



5/2 



1 + 7 



1751 7 
l336~ ~ 1 U 



(2.7) 



Using the post-Newtonian law (2.4) truncated at 1PN order we further deduce 

96 Gm Klo I" / 2591 



1/7 
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1 + 7 
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(2.8) 



These expressions can be obtained alternatively using the center-of-mass 3.5PN equations 
of motion for general orbits [331 S3]- By substituting them into Eq. (2.2) we obtain the 
complete acceleration for quasi-circular orbits at 3.5PN order as 



dv 

dt 



-co 2 X 



32 Gm 



5/2 



1 + 7 



743 11 



336 



-v 
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(2.9) 



where the radiation reaction force up to 3.5PN order is explicitly exhibited. As for the 
velocity it is given with the same accuracy by 



v 



ruX — 



64 Gm 



1/7 



5/2 



1 + 7 



1751 7 
"336" ~ 1 P 



n 



(2.10) 



Those are the coordinates for which the metric is asymptotically Minkowskian far from the source. 
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FIG. 1. Definition of angles and polarization vectors. AT" points towards the observer, whose 
spherical coordinates are (0, <£), with B = ir/2 corresponding to the orbital plane of the binary; O is 
also the inclination of the orbit on the sky from the observer's viewpoint. We define then P = — e$ 
and Q = bq. The (n, A, £) triad is displayed, <p being the orbital phase, and 4> = f — $ + 7r /2 
being the orbital phase relative to the direction of the ascending node of the orbit, which is also 
the direction of P. 



III. GRAVITATIONAL WAVEFORM OF NON-SPINNING BINARIES 



A. General definitions 



The gravitational waveform propagating in the asymptotic regions of an isolated source 
is defined in a radiative (Bondi-type) coordinate system = (cT, X) by [35] 

AG + °° if 2£ 1 

h 7 = ^ViMN)J2^\ N ^ U kiL-2(T R ) --^—^N aL ^e ab[k V l)bL - 2 (T R )j , (3.1) 

where terms of order 0(1/ R 2 ) are neglected. Here R = \X\ and N = X j R are the distance 
and the direction of the source. The Minkowskian retarded time is denoted Tr = T — R/c. 
The transverse-tracefree (TT) projection operator reads Vijki = Vi(kPi)j — \PijVu where 
Vij = Sij — NiNj is the projector orthogonal to the unit direction N. The waveform (3.1) 



is parametrized by two sets of radiative symmetric-trace-free (STF) multipole moments, 
of mass-type, Ul, and of current-type, Vl, which are functions of the retarded time Tr. 
The notation for STF tensors and multi-indices such as L — i\ • ■ ■ if (where £ denotes the 
multipole order) is the same as in Ref. [20J (Paper I) where it is fully detailed. 

We now introduce two unit polarisation vectors P and Q forming with N an orthonormal 
triad, as indicated on Fig. [1} Representing the observers's direction in spherical coordinates 
(0,$), with O = 7r/2 corresponding to the orbital plane of the binary, our choice for the 
polarization vectors (which is not unique^) is P = — e$ and Q = e e , where e e = dN/dQ 



5 Beware that other authors use different conventions for the polarization triad. In Ref. |19j . this results in 
an overall sign difference with us for both /i + and h x . 
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and e$ = 9 A/"/ (sin 0<9$) are the standard transverse basis vectors of spherical coordinates. 
With respect to the polarization vectors P and Q, the two polarization states are 



h. 



\ {PiPi - QiQj) hj? 
liPiQj + PjQjhJ? 



(3.2a) 
(3.2b) 



We shall henceforth pose h = h + —ih x . For quasi- circular orbits the complex wave-amplitude 
h will be a function of the orbital phase if = judt and of the spherical angles (9, $) 
associated with N. Next, we apply a decomposition in terms of spin-weighted spherical 
harmonics Y_^(Q, $)|^] This yields the definition of the spherical modes h £m of the waveform 

as 

+oo e 

h = J2Y, h £m Y^{®^). (3.3) 

1=2 m=-l 

From the orthogonality properties of the Yf™'s, the separate modes h lm are extracted from 
h by the surface integral (with solid angle element d 2 Q = sin 0dOd$) 

h im = J tfQhY^Q, $) , (3.4) 
where the overbar means the complex conjugation. We can thus insert the waveform de- 



composition (3.1) in terms of STF radiative moments and perform the angular integration, 
to obtain 
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(3.5) 



where U £m and V £m are the radiative mass and current moments in standard (non-STF) 
form [19j H5], related to the STF radiative moments by 



U 



4 l(£+l)(£ + 2) . 

a L U L 



2£(£-l) 



V 



£{e + 2) 



£\ V %t + !)(£- 1) 



_ im t t 

Oi T Vi 
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(3.6a) 
(3.6b) 



Here we have introduced the tensor a 1 ™ defined by the decomposition of the STF product 
of £ unit vectors N, i.e. Nl in our notation, into ordinary spherical harmonics, namely 



Nr 



£ ai m r £m (e,$) 



(3.7) 



See Ref. |19| ; our convention for the spin- weighted spherical harmonics is the same as in Paper I. 
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B. Mode separation for planar binaries 



Let us now show that for planar binaries, the mode h is in fact entirely given by 



the "mass" contribution [first term in (3.5)] when £ + m is even, and by the "current" 
contribution [second term in ( |3.5[ )] when t + m is odd. Notice that this is valid in general 
for non-spinning binaries, regardless of the orbit being quasi- circular or elliptical, or for 
spinning binaries with spins aligned or anti-aligned with the orbital angular momentum. 
The important point is only that the black holes must be non-precessing so that the motion 
is planar. This mode separation has been already used in Ref . [19] , but we explain here how 
to deduce it from an argument of parity invar iance. If an observer with space-coordinates X 
in a spatial grid {X 1 } looks at the gravitational- wave signal at a given time T, he finds that 
the source particles lie at some positions y A and possess certain velocities v A , for A = 1,2. 
Now, a second observer located at X' = —X and looking at a parity-reversed system with 
positions y' A = —y A and velocities v' A = —v A must observe the same waveform hj?, by 
parity invariance of the Einstein field equations. This translates intoj^] 



hJ?(-X, T; -y A , -v A ) = hJ?(X, T; y A , v A ) . (3.8) 



Because the multipolar decomposition (3.1) is unique, this result allows us to recover the 
basic built-in parity properties of the mass and current moments Ul and Vl, regarded as 
functionals of the matter variables y A , v A : 

U L (T; -y A , -v A ) = (-1)" U L (T; y A , v A ) , (3.9a) 
V L (T;-y A ,-v A ) = (-l) i+1 V L (T;y A ,v A ). (3.9b) 



To see how Eq. (3.8) reflects on the two polarizations, we remark that the polarization vector 
P = N x £/\N x £\ transforms as P' = —P, while Q = N x P is clearly invariant, due 
to the fact that £ is itself invariant and N' = —N. Those considerations applied to the 



complex polarization h = h + — ih x , with the definition (3.2) for h + and h x , show that, in 
the center-of-mass frame 

h(-X,T]-x,-v) = h(X,T;x,v), (3.10) 

where x = y\ — t/2 and v = dx/dt. 

For planar motions, it is particularly convenient to rewrite the previous relations in the 
spherical coordinates associated with the frame {X 2 }. The triad {N,P, Q} is specified by 
the inclination G and azimuthal angle $ of the unit vector N. Similarly, the mass-point 
location is determined by the pair of angles (6 = tt/2, (p), where ip is the orbital phase. The 
counterparts of those quantities for the reversed observer are O' = 7r — 0, $' = $ + 7r, 
9' = 7r/2 and ip' — (p + 7T. Note that the components of the relative position x and velocity 
v are solely functions of the azimuthal phase (p on the trajectory r = r(<p). Now as is 
clear from the geometry of the planar binary, the invariance under rotation implies that 
h(R, T, 0, $, (p) depends on $ and <p> only through their difference ip — Hence we define 
as a convenient phase variable <fi = (p — $ + 7r/2, which is actually the orbital phase whose 
origin is chosen to be the direction of the ascending node P (see Fig. [T]). We shall regard h 
as a function of and 0, denoted by h(Q,4>) hitherto, ignoring dependences in R, T since 



7 In the non-precessing case, the aligned or anti-aligned spins are only characterized by their magnitude, 
which is a constant parameter that does not play a role here. 
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they do not play a role here. Putting those results together, we arrive at the important 
identity 

h(-K-e,(j)) = h(e,(f)). (3.11) 

Inserting the definition (3.4) of h into both sides of the latter equation and noticing [4*6] 

that Y e ™,(ir - 0,$ + 7r) = (-l) m+ ^Ff 2 " m (6, $), we conclude from the uniqueness of the 
decomposition in spin-weighted spherical harmonics that 

h tm (<p) = (-iy +m h e '~ m (ip + it) . (3.12) 

The known dependence of the mode on the azimuthal angle, h lm {(p) oc e~ imip , then shows 
that 

h tm (tp) = (-l) e lf'~ m (tp). (3.13) 

In the end, we substitute for h im its expression ( 3.5[ ) in the equality above. The tensor a£ m 
introduced in Eqs. (3.6) satisfies a e l~ m = (— l) m a^ m J^] which implies that the multipolar 

coefficients U tm and V tm satisfy themselves the relations U ' m = (— l) m U lm and V ' m = 
(— l) m y' ra , Finally, we find that these multipole coefficients must be such that U lm = 
(— \y+ m jj em anc i v £m = (—iy+ m + 1 v e ' m , hence we have necessarily U im = when £ + m 
is odd, and V £m = when i + m is even. We therefore conclude that there is no mixture 
between the mass-type and current-type contributions to the modes, in the sense that 

h lm = U lm when i + m is even , (3.14a) 

h lm = i V tm when i + m is odd . (3.14b) 

Let us stress again that this result holds only when the orbit stays planar and, in particular, 
does not apply to precessing binaries. However, it is valid for quasi-circular as well as 
elliptical orbits. 



IV. EXPRESSION OF THE DOMINANT QUADRUPOLE MODE 

In the following sections, we shall compute the dominant quadrupole mode (£, m) = (2, 2) 
at 3.5PN order. From the results (3.14), the mode h 22 depends only on the radiative mass- 
type quadrupole moment ?7yF] 



22 _ rr22 _ V&G 2 '2 

v'2/.V Jfc* ' " * 



U 22 = ~^» 22 U l3 . (4.1) 



Following the multipolar-post-Minkowskian (MPM) formalism [251429] . we shall express the 
radiative quadrupole moment in terms successively of canonical moments {Ml, Sl} and 
then source moments {1^, Jl,Wl, X L ,Y L , Z L } for general sources at 3.5PN order. Among 



Sec Eq. (2.12) in [45]. The notation used in Refs. QHH] is related to ours by y e L m = (2 ^, )!! af 
The STF tensor otfj therein is explicitly given by 

with the brackets surrounding indices relerring to the STF projection. 
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the source moments, {II, Jl} are the most important ones because they arise at Newto- 
nian order, while {Wl, X l , Y l , Z l }, which are called the gauge moments, induce corrections 
starting at high 2.5PN order. 

Using the "selection rules" developed in Paper I, we can identify all types of terms to be 
computed up to quadratic order (which is sufficient for our purpose). Let £ max [Aj, Bk] be the 
maximum multipolar order beyond which Ul cannot contain products of the two multipole 
moments Aj and Bk or their time derivatives (or possibly time anti-derivatives), at the 
3.5PN level. Here Aj and Bk can be any moments chosen among the set of source moments 
{II, Jl, Wl, X l , Y l , Z l }. Thus when £ > £ max [Aj, B K ] any term in Ul constructed with the 
product of moments Aj and Bk has a post-Newtonian order strictly higher than 3.5PN and 
is negligible here. On the other hand when £ max [Aj, Bk] < 2, the corresponding interaction 
is not allowed to enter the quadrupole moment Uij which has £ = 2. When £ raax [A J , B K ] ^ 2, 
it is equal to the number of (free or dummy) indices to be shared between: (i) Aj and Bk 



k , or n 



A a 



J-ij 



BhK 



if both moments are of the same parity, i.e. j + k = 

and Eiab (with contractions on indices a and b) if both moments are of different parity, 



:[Aj,B 



bK-1 

i.e. 



j + k 



■\Aj, Bk] + 1. The values of £ max [Aj, Bk] larger than or equal to 2 are listed in 



Table [I] for the various possible choices of Aj and Bk- Take for instance the case of Ij and 
Ik- By dimensionality we can add the powers of 1/c and find that this term is necessarily 
of the type (p and q being the number of time derivatives) 



Ui 



G> 



_^ T (P) ril) 

- k +3-£ -J A K 



(4.2) 



The contribution of this term in the waveform (3.1) is easily seen to be of post-Newtonian 
order 0(l/c> +k+1 ) and hence we require j + k = 6 for a 3.5PN term. On the other hand, 
we have said that £ max = j + k in this case, hence we have £ max = 6 in Table [j] for a 



contribution like (4.2). Similar selection rules can be made for the radiative current moment 
Vl but here we need only to compute the mass quadrupole moment. From Table |TJ it is 
straightforward to obtain the set of quadratic interactions contributing to U^ within our 
accuracy. No new cubic interaction (between three source moments), besides the already 
known 3PN tail-of-tail term, enters the mass quadrupole Uu at the 3.5PN order. 



Aj 


Bk 


Ik 


Jk 


W K 


Xk 


Y K 


Zk 


Ij 


6 


4 


4 


2 


4 


2 


Jj 


4 


4 


2 




2 


2 


Wj 


4 


2 


2 




2 




Xj 


2 












Yj 


4 


2 


2 




2 




Zj 


2 


2 











TABLE I. Values of £m&x[Aj, Bk] at the 3.5PN order for the various possible choices of Aj and 
Bk- Because factors may always be reordered in products, we have l max [Aj,BK] = £m&x\Bj-, Art] 
and the table is trivially symmetric. 



The long calculation following the MPM algorithm [25] (in the slightly modified version 
of Ref. [27]). is implemented in a systematic way by means of the tensor package xTensor 
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developed for Mathematica [47], and yields various types of terms: 



U, 



mem 

ij 



(4.3) 



In a first stage we express each of the different pieces by means of the canonical moments 
{Ml, Sl}. The instantaneous piece up to 3.5PN order is 
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(4) 
3) be 
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(4.4) 



where the angle brackets ( ) around indices % and j denote their STF projection (for conve- 
nience we do not indicate the neglected remainder 0(l/c s )). The hereditary tail integrals 
involve the dominant tail term at 1.5PN order [26] and the tail-of-tail term at 3PN order 
EH1: 



t rtail 
ij 



2GM 



dr 
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H In 
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Tr-t 
2r 



124627 
44100 



M^(r), (4.5) 



where To denotes an arbitrary constant which will enter in fine into some unphysical fre- 
quency scale (6.2) in the waveform. The memory-type hereditary integrals, which constitute 
the third piece in Eq. (4.3), contribute here at orders 2.5PN and 3.5PN. Our expression ex- 

ESH5D], 



tends the classic works on the dominant 2.5PN non-linear memory effect j26J E 
and is in agreement with the recent computation of the non-linear memory up to any post 
Newtonian order in Refs. [5T| 521. We have 
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• (4-6) 



In a second stage of the general formalism, we must express the canonical moments {Ml, Sl} 
in terms of the six types of source moments {II, Jl, Wl, Xl, Yl, Zl}- For the present com- 
putation we need only to relate the canonical quadrupole moment My to the corresponding 
source quadrupole moment Jy up to 3.5PN order. We obtain 
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(4.7) 



The current quadrupole canonical moment and mass octupole canonical moment My*, 
agree with the corresponding source moments up to 2PN order, and for this computation it 
is sufficient to replace these by their source counterparts. 



V. SOURCE MULTIPOLE MOMENTS OF COMPACT BINARIES 



We now feed the previous expression of the radiative quadrupole moment t/y with explicit 
results for the source moments in the case of non-spinning compact binaries on quasi-circular 
orbits. The general formulas for the six sets of source moments {II, Jl, Wl, Xl, Yl, Zl} as 
explicit integrals over the matter distribution in the source and the gravitational field were 
obtained in Ref. [22] • The mass-type source moments expressed in terms of elementary 
retarded potentials up to 3.5PN order are given by Eqs. (3.3)-(3.6) of Ref. [31] • Reducing 
the latter source quadrupole moment ly for compact binaries, in the frame of the center-of- 
mass and for quasi-circular orbits, yields 
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where x l is the orbital separation and v % = dx l /dt is the relative velocity. Resorting also to 

defined by Eq. (2.3), the "conservative" coefficients A and 



the post-Newtonian parameter 7 
B in (5.1 ) read 
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(5.2b) 



The results at 3PN order were already known [30, 31 J . Note the appearance at 3PN order of 
two arbitrary scales: the first one is given by r$ = ctq, where tq is the same scale as in the 



tail integrals (4.5); the second one is a Hadamard regularization scale r' due to the presence 
of the point mass singularities. The scale r is associated with an irrelevant choice of the 
origin of time in the radiative zone, whereas the scale r' will disappear from our final results, 



being cancelled by the corresponding r' in the equations of motion; see Eq. (2.4). The new 
correction term at 3.5PN order only affects, for circular orbits, the "radiation reaction" 
coefficient Crr which is found to be 
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(5.3) 



We give the result for the complete 3.5PN mass quadrupole moment in the case of generic 
orbits reduced to the center-of-mass frame in Appendix [Aj 

In addition to the 3.5PN mass quadrupole iy we will also need the expression of the mass 
monopole or ADM mass M (notice that / = M) up to 2PN order since it is required in the 



computation of the tail terms (|4.5|). We have, for circular orbits: 

V V , x c, „ / 1 



M = m 



7 + o(7-^)7 2 



+ o 



(5.4) 



Furthermore the current dipole moment Jj is required up to 1PN order, since it appears in 



the 2.5PN term of Eq. (4.4) (recall that Si = Ji in our approximation): 
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The other source moments II and Jl are required at Newtonian order only, and at that 
order we have the following general formulas, valid for all multipolar indices I ^ 2: 
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and X 2 = ^ = ^(1 - A) where A = = ±y/l - 4z/; see also Footnote 24 of Paper I 

for an alternative expression of the function se{v). 

Among the required gauge-type source moments {Wl, Xl, Yl, Zl}, the only one which is 
to be controlled at 1PN order is the monopole W, since it enters the 2.5PN term of Eq. (4.7). 



However it turns out to be zero for quasi-circular orbits. The other needed moments are 
merely Newtonian. Their full list is as follows: 
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Finally, when computing the waveform, we need the multiple time derivatives of all these 
moments; they are computed by order reduction using the acceleration for quasi-circular 
orbits given at 3.5PN order by Eq. (2.9). Furthermore, we need to perform many scalar 
products between the velocity and the polarization vectors (N, P,Q); these are computed 
using the 3.5PN expression of the velocity as given by Eq. ( |2.10 ). Note also that the (tails 
and non-linear memory) hereditary terms given by Eqs. (4.5 ) and (4.6 ) are treated by exactly 
the same method as in Paper I. 



(-3 + 16v)x {i x j) -?mi/(l- 3i/)rW> + O 



VI. DOMINANT QUADRUPOLAR MODE AT 3.5PN ORDER 

We shall now perform a change of phase variable, from the actual orbital phase related 
to the orbital frequency by <p = J udt, to a new phase variable ifi, defined in such a way 
that it absorbs most of the logarithmic contributions in the frequency coming from the tail 
terms. Extending Refs. [HH \T7\ 120] . we find that the new phase variable takes remarkably 
the same expression as the one used in Paper I: 



IGMuj /u\ 

c 3 \Uo) 



^ = V — ln( — ). (6.1) 



The 2PN corrections in the ADM mass M, as given by Eq. (5.4), are exactly the ones 



needed for the logarithmic cancellations at our extended 3.5PN order. Here ujq is a constant 



frequency scale related to the constant r appearing in the tail integrals (4.5) by 



11 -7E 



UJ = ^— , (6.2) 
4r 

with 7e = 0.577- • • denoting the Euler constant. It is known [TU] that this frequency scale 
is irrelevant: a rescaling of Uq is indeed equivalent to a shift of the binary's instant of 
coalescence, which can be absorbed into a redefinition of the origin of time in the radiation 



zone. Physically, the logarithmic term in Eq. (6.1) expresses the fact that the tails induce 



a small delay in the arrival time of gravitational waves. When viewed as a modulation of 



the orbital phase as in (6.1), this term represents in fact a very small correction to this 
phase, being of order 4PN when compared to the dominant phase evolution which is at the 
"inverse" of 2.5PN order (i.e. at order c +5 ), and is negligible with the present accuracy. 
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Replacing in terms of the new phase variable ip, we perform the appropriate Taylor 



expansion of the waveform and apply the mode decomposition (3.3). We also use the fact 
that, from the planar geometry of the problem, h = h + — ih x depends only on the inclination 
angle G and the difference ip — <£> between the "absolute" phase ip and the azimuth $ of the 



vector N (see also Sec. Ill and Fig. [TJ. For instance, like in Eq. (3.11) but with ip playing 



now the role of <p, we can introduce the phase <p = ip — $ + 7r/2 defined relatively to the 
direction of the ascending node P. Choosing the latter phase as integration variable instead 
of $, and using the known dependence of the spherical harmonics on the angle 0, we obtain 

h tm {ip) = (-i) m e~ imip [ an h{e, $) ri 2 m (e, $) , (6.3) 



exhibiting the azimuthal factor e appropriate for each mode. Alternatively, we may 
access directly a given h mode, once the requisite radiative multipole Ul or Vl known, by 



using the contraction formulas (3.6), which translate into (4.1) for h 22 . Finally, to present 
our result we pose 

h tm = 2G ™" 2 X yJ^H^e-W , (6.4) 
and express it entirely with the help of the gauge invariant post-Newtonian parameter x 



related to the orbital frequency by Eq. (2.5). We find for the dominant (2,2) mode up to 



3.5PN order (extending the 3PN result of Paper I for this mode): 
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We have verified that the above expression is in full agreement in the test mass limit where 
v — > 0, with the result of black-hole perturbation theory as reported in the Appendix B of 
Ref. [23] or, in a form more suitable for comparison, in Eq. (4.9a) of Ref. [24] . 

As it is the dominant quadrupole (2, 2) mode that is determined with the best precision 



by numerical relativity, the 3.5PN result (6.5) should become important for high-accuracy 



comparisons between numerical and analytical waveforms. See Paper I for all the other 
modes that contribute up to 3PN order. The completion of the remaining modes at 3.5PN 
order presents some new difficulties; in particular the (2, 1) mode necessitates the difficult 
extension of the current quadrupole moment Jij up to 3PN order. This is left for future 
work. 
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Appendix A: 3.5PN mass quadrupole moment 



We give our result for the complete source mass quadrupole 1^ up to 3.5PN order, for 
general non-circular orbits, but reduced to the center-of-mass frame. Defining: 
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Those expressions were already given at 3PN order in Eqs. (5.10) of [21] • The constants r 
and r' match the ones introduced in Eqs. (5.2). The parts v4rr, and Crr corresponding 
to the radiation-reaction contributions and containing the new 3.5PN terms are: 
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Note that we defined -4rr, -Brr, and Crr as being dimensionless, but that ^4rr and Crr 
start at 2.5PN order, while £>rr only starts at 3.5PN order. The expressions (A2)-(A3) 
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reduce to Eqs. (5.2)-(5.3) in the case of quasi-circular orbits. Notice also that we replaced 
the Hadamard regularization ambiguities £, ( and k in the 3PN terms by their values, which 
can be found for instance in Ref. 
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